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Abstract 

The Legendre flow equation, a version of exact Wilsonian renormalization group (WRG) equa- 
tion, is employed to consider the power counting issues in Nuclear Efl'ective Field Theory. A WRG 
approach is an ideal framework because it is nonperturbative and does not require any prescribed 
power counting rule. The power counting is determined systematically from the scaling dimensions 
of the operators at the nontrivial fixed point. The phase structure is emphasized and the inverse 
of the scattering length, which is identified as a relevant coupling, is shown to play a role of the 
order parameter. The relations to the work done by Birse, McGovern, and Richardson and to the 
Kaplan-Savage- Wise scheme are explained. 
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I. INTRODUCTION 



Conventional nuclear theory is based on force potentials. There have been much progress 
and we now have very accurate ones such as NijmegenQ], Argonne VisPl! and CD-Bonn 3,01 
potentials which fit nicely to about 3,000 nucleon-nucleon scattering data with energies 
up to 350 MeV. The potential models have matured to the demand of precise numerical 
calculations. 

Potential models however have serious drawbacks. Even though they are very precise, 
after all, they are just (semi-phenomenological) models. It is not obvious how these models 
are related to QCD, the fundamental theory of strong interactions, nor how to improve 
them in a systematic manner. Proper treatment of inelastic processes seems formidable. In 
addition, the potential approach is known to suffer from the so-called off-shell ambiguities, 
which may cause serious problems for a correct understanding of multi-nucleon systems. 

Nuclear effective field theory (NEFT)P, 0] is expected to _be a promising^ alternative. It 
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and linked to 



is based on very general principles of quantum field theory 
QCD through chiral symmetry, while, instead of quarks and gluons, only the relevant degrees 
of freedom (nucleons, and optionall y p ions, deltas, etc.) at low energies are considered. 



Chiral perturbation theorvf yPT) |l3l. Il4| . a prominent example of EFT in which only meson 



degrees of freedom are treated explicitly, has been applied very successfully. In short, NEFT 
aims at a similar success. 

Since a Lagrangian of EFT contains an infinite number of operators allowed by the sym- 
metry, one needs a power counting rule to give orderings to an infinite number of Feynman 
diagrams generated by the Lagrangian. Once it is specified, physical quantities can be cal- 
culated in a systematic expansion in powers of Q/Aq, where Aq is the cutoff, above which 
the EFT is not valid, and Q (< Aq) is the typical momentum of the process one is inter- 
ested in. In the case of perturbative systems such as meson-meson scattering, the power 
counting which is determined on the basis of naive dimensional analysis(NDA) 1^ is known 
to give correct results. NEFT is, on the other hand, inherently nonperturbative. (There 
are bound states, the nuclei!) In addition, it is a fine-tuned system; the inverses of the 
nucleon-nucleon scattering lengths l/a'-^'^"-' ~ —8 MeV, and ~ 40 MeV are much 

smaller than Aq ~ ~ 140 MeV for the pionless NEFT. (If we consider the pionful NEFT, 
Aq should be larger, presumably Aq ~ 500 MeV, and the system should be considered more 
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fine-tuned.) Moreover deuteron is a shallow bound state with binding energy of 2.2 MeV. 
We need a power counting scheme which takes into account the "unnaturalness." 

In this paper, we consider a systematic way of determining the power counting for such 
nonperturbative and fine-tuned systems on the basis of Wilsonian RG analysis. As we will 
explain below, our approach is very general and not restricted only to NEFT. 

At present, there are basically two power counting schemes proposed for NEFT. In Wein- 
berg's scheme js,^], one identifies the sum of irreducible diagrams in the expansion based 
on NDA as the "effective potential," and solves the Lippmann-Schwinger (L-S) equation 
with it. The "effective potential" is iterated to all order, i.e., all operators are treated non- 
perturbatively. Although conceptually attractive and simple in numerical implementation, 
Weinberg's power counting has been shown to be inconsistent: the divergences that arise 
in the leading order calculation cannot be absorbed by the leading order operators in cer- 



tain partial waves 
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17| . The other power counting scheme proposed by Kaplan, Savage 



and Wise (KSW) [JJ, be gins with a renormalization group (RG) behavior at a nontrivial 
fixed point, which corresponds to a system with infinite scattering length. In contrast to 
Weinberg's scheme, only the nonderivative contact operator is treated nonperturbatively 
while others are treated perturbatively. Although the KSW power counting is a consistent 
scheme in the sense that results are independent of subtraction point at every order of the 
expansion, serious problems have been reported. Fleming, Mehen, and Stewart fFMS)[2ol| 
calculated the nucleon-nucleon scattering amplitudes up to next-to-next-to leading order in 
a NEFT with pions included ("pionful" NEFT), implementing the KSW power counting, 
and found that the scattering amplitude in the ^Si-'^Di channel shows no con verg ence at 
all. A possible explanation for the success and the failure has been given in Ref . |2ll . |22 1 . 

The present status may be summarized as that NEFT works "pretty well if one follows 
a patchwork of power counting rules j^]." Actually, in spite of the problems mentioned 
above, many of higher order calculations in both power counting schemes have been found 
to be successful at least numerically. It might be however too early to say that the problems 
are settled. One should keep in mind that, besides the numerical agreement with data, a 
thorough understanding of power counting in a systematic and consistent way is extremely 
important. It is this issue that we address in this paper. 



We believe that Wilsonian Renormalization Group (WRG) 



24 |25| approach can provide 



the key. This approach allows us to analyze the behavior of the operators (or the corre- 
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spending coupling constants) by assuming only the symmetry and the relevant degrees of 
freedom. Most importantly to our present purpose, WRG analysis is nonperturbative and 
does not require any prescribed power counting scheme. 

Our basic assertion is that it is anomalous dimensions that determine the power counting. 
We are going to show in this paper that it is very natural to determine the power counting 
of an operator from the corresponding scaling dimension (the sum of the canonical and 
the anomalous dimensions). Our reasoning is based on the very general idea that the power 
counting should be determined by dimensional analysis. Furthermore, our WRG approach is 
applicable both to trivial and nontrivial fixed points; near the trivial fixed point, it coincides 
with NDA, while it explains why the KSW power counting is successful to some extent near 
the nontrivial fixed point. The real two-nucleon system happens to be close to the nontrivial 
fixed point, but the WRG approach has a wider scope. 

Wilsonian approach has been applied by Birse, McGovern, and Richardsonj^] to nonrel- 
ativistic two-body scattering to obtain a power counting rule. They considered a Wilsonian 
RG equation for an effective potential, and identified two fixed points; a trivial one and a 
nontrivial one. Then they determined how to organize the terms in the effective potential by 
considering the linearized RG equation around the fixed point. Although the potential they 
considered is general, the physical meaning of energy-dependence of it in the Schrddinger 
equation is obscure. Furthermore it does not seem to be very clear how the symmetry of 
the theory is implemented in the potentials. Since this potential approach is similar to tra- 
ditional potential model approach, it might share some of their problems. It is desirable to 
understand the power counting in a completely field theoretical framework to benefit from 
all the good features of EFT. The relation between their work and the present work will be 
explained in Sec.|V|in detail. 

A WRG analysis in field theory requires all possible operators consistent with the as- 
sumed symmetry. In the previous paper j^l, we emphasized the role played by "redundant 
operators," which may be eliminated by field redefinitions. We then explicitly showed that 
the nonperturbative off-shell amplitude can be renormalized only by properly treating re- 
dundant operators in the WRG analysis. However, the method of the calculation is rather 
straightforward and seems applicable only for simple cases. A more powerful method is 
needed for general cases. 

In this paper we study the RG flow for the pionless NEFT employing Legendre flow 
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equation, which is one of the implementations of Wilsonian RG. We obtain the RG flows for 
^5*0 and ^Si-'^Di channels to order 0{p^). We identify two phases in our RG flow for both 
channels; the strong coupling phase and the weak coupling phase. We also show that the 
nontrivial flxed point is on the phase boundary of our RG flow. The scaling dimension of 
the operators around the flxed points are calculated. We flnd that the anomalous dimension 
of the operators around the nontrivial flxed point is large. We argue that the scattering 
length is identifled with the order parameter which characterizes these phases. 

The pionless NEFT is an interesting theory as an application of nonperturbative RG 
equation, which is, in most cases, very difficult to treat. Thanks to the nonrelativistic 
feature of the theory, it allows us to solve (approximately) with a simple truncation of 
the space of operators, though the justification of the truncation only comes from actually 
enlarging the space. This theory also provides a simple example of nontrivial fixed point 
and the existence of a bound state. We hope that analyzing this theory gives some insight 
into the use of nonperturbative RG to get the information about bound states. 

This paper is organized as follows: in Sec. |nj we explain why the scaling dimensions are 
important in determining power counting. The Legendre fiow equation for nonrelativistic 
systems is introduced in Sec lIIII We present the detailed analysis of pionless NEFT and its 
physical implications in Sec. IIVI The relation to the work done by Birse et al. is clarified 
and comments on the other schemes, especially on the KSW scheme, are made in Sec. 
Finally in Sec. IVH we summarize our study, and discuss future prospects. Appendix |X] 
provides some technical information about the cutoff function used in the Legendre flow 
equation. In Appendix |Bl we derive the RG equations from the two-nucleon scattering 
amplitudes. 

II. POWER COUNTING AND RENORMALIZATION GROUP 

In this section, we consider what power counting is and its relation to renormalization 
group in general terms. 

The most basic idea behind power counting is the order of magnitude estimate based 
on dimensional analysis. The order of magnitude of a physical quantity may be estimated 
by an appropriate combination of dimensionful parameters of typical scales. The period T 
of a simple pendulum with length L and mass M in a uniform gravitational fleld with the 
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gravitational acceleration g is estimated as T ~ ^jL/g. The dimensionless coefficient is 
expected to be of order one, the idea of "naturalness." (In the pendulum example, it is a 
bit large, 27r.) 

There is only one important scale in a simplest (relativistic) field theory; the physical 
cutoff Aq. The mass m is related to Aq by a very large dimensionless constant ^ as m = Aq/^. 
(This is a necessary fine-tuning to have a sensible quantum field theory.) Any quantities 
may be measured in units of Aq. The interaction Lagrangian in spacetime dimension D may 
be written as an infinite sum of the interactions of the form 

^^nt = -J2^0.{x), (2.1) 

i 

where di is the (canonical) dimension of the operator Oi{x) and g^i is the (bare) dimensionless 
coupling constant. (It is sometimes useful to say that the (dimensionful) coupling goi = 
(7oi/Ao'~'^ has dimension D — di.) Classically (and near the trivial fixed point, see below) 
this dependence on Aq determines how large the contribution of Oi is. One may expand the 
contribution in powers of (Q/Aq), where Q is the typical energy /momentum scale of interest, 
according to the canonical dimensions of the operators. This is a common situation, and, 
with a twist of chiral symmetry, true for xPT^. The so-called NDA applies to such cases. In 
some cases, however, quantum fluctuations change the situation drastically. 

In a WRG analysis, the cutoff is lowered to A < Aq by integrating out the fluctuations 
with momentum A < p < Aq, and thus the interaction Lagrangian may be replaced by 

^- = -E|^^^^(^)- (2-2) 

i 

The coupling gi{A) now depends on the "floating" cutoff A so that the physical quantities 
do not depend on A. The behavior of gi{A) contains information about how the quantum 
fluctuation modifies the importance of the operator O^. The differential equations, 

t^ftte). *...(^). (2.3) 

^ In xPT, the masses of the Nambu-Goldstone bosons play a special role as symmetry breaking parameters, 
which arc not related to the cutoff Aq. It is not trivial how to treat them. The standard way is to treat 
~ Q^, which is necessary to be consistent with the meson pole structure of the amplitudes jl^. In the 
pionless NEFT. which wc consider in this paper, there is no such a symmetry breaking parameter (we 
ignore isospin breaking), and thus no subtlety associated to it. 
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are called RG equation, and determine the RG flow in the space of all coupling constants. 

An important concept in the RG analysis is a fixed point of the RG fiow, which is a 
solution of Pi{g*) = for all i, where the coupling constants stop running. Suppose we are 
interested in the behavior of a theory close to a fixed point. The behavior of the RG flow 
near a fixed point can be examined by considering the linearized RG equation obtained by 
substituting gi = g* + 6gi into the RG equation keeping only the linear terms in 6gi, 

^^5g, = A,{g*)5g„ (2.4) 

where Aij{g*) = {d[3i/dgj)\g*. By diagonalizing Aij{g*), we have 

dui 



dt 



UiUi, (2.5) 



where Vi is the eigenvalue and Ui is the corresponding eigenvector, which may be immediately 
integrated as, 

u,(A)=^,(Ao)(^) \ (2.6) 

The exponent of Aq, i^j, is called scaling dimension of the coupling. Depending on the sign 
of the scaling dimension, interactions are divided into two groups. An operator with its 
coupling having z/j < is called irrelevant. On the other hand, that with z/j > is called 
relevant. If Vi = 0, the operator is called marginal. The scaling dimension is the quantum 
counterpart of the (canonical) dimension. At the trivial fixed point, g* = for all i, it agrees 
with the canonical one, Ui = D — di, i.e., 

while at a nontrivial fixed point it may be quite different. It should be emphasized that the 
scaling dimensions are not prescribed but determined by the theory (and the fixed point) 
itself. This is the beauty of the WRG approach. 

Let us consider what happens near a nontrivial fixed point more in detail. Suppose that 
coupling constant gi{A) is written as 

where Cik is a small constant. First of all, it is important to note that the theory is scale 
invariant at the fixed point. It means that the theory has no typical scale, though it starts 



with the physical cutoff Aq. Actually, the coupling gi{A)/ A^'~^ ~ g*/A'^^~^ is independent 
of Aq. Therefore, the value of the fixed point itself does not contribute to the power counting, 
because, as we emphasized above, power counting is based on dimensional analysis, but a 
scale invariant theory does not have a scale! It is the Aq dependent part (the second term 
in (j2.8|) ) that contribute to the power counting. 



From the dependence on Aq, one sees that the k-th term in the sum behaves like a coupling 
with dimension z/^. 

In the literature, many authors seem to think that the lowest order contact operator 
in NEFT should be treated "nonperturbatively," because the value of the corresponding 
coupling constant is large at the nontrivial fixed point (the first term of (|2.9|) ). As we will 
show in Sec. IIVI however, it should be treated "nonperturbatively" because it is relevant at 
the nontrivial fixed point. We believe that this difference is crucial in understanding the 
power counting for other operators. 

Let us summarize what we explained in this section; power counting is in essence the 
order of magnitude estimate based on dimensional analysis. But the quantum fluctuations 
change the dimension drastically in some cases. A WRG analysis takes into account the 
quantum fluctuations in a controlled way, and gives rise to a RG flow, which is characterized 
by fixed points. Especially near a nontrivial fixed point, the scaling dimension may be quite 
different from the canonical one. It is therefore natural to consider the power counting based 
on the quantum scaling dimensions. The important point is that it is not the value of the 
nontrivial fixed point, but the scaling dimensions that contribute to the power counting. 

Finally, we would like to make a comment on the use of momentum cutoff in EFT. 
Although it might not be widely understood, the momentum cutoff scheme does not break 
the EFT expansion. (Of course it is much simpler with dimensional regularization.) The 
loop contributions which (partially) cancel the Aq in the denominator may be absorbed in 
the renormalized couplings of lower order. Thus, when expressed in terms of renormalized 
couplings, EFT with momentum cutoff has a sensible expansion as that with dimensional 
regularization. 




(2.9) 
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III. LEGENDRE FLOW EQUATION FOR NON-RELATIVISTIC THEORY 



Legendre flow equation 
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which we use in our analysis, is one of the implemen- 
tations of WRG equation. It is formulated as a RG equation for the infrared (IR) cutoff 
effective action called effective average action. In contrast to the conventional effective ac- 
tion in which all the fluctuations are taken into account, the effective average action includes 
only the quantum fluctuations with momenta larger than some cutoff A. An infinitesimal 
change of an effective average action with respect to an infinitesimal change in A is expressed 
in a differential equation. Starting from an arbitrary bare action at ultraviolet (UV) scale 
Aq, the effective averaged action successively includes the lower momentum fluctuations as 
the cutoff is lowered, approaches the conventional effective action in the limit of A = 0. 

In the case of relativistic theory, it is necessary to define the flow equation in Euclidean 
space because a cutoff must be imposed on all four components of the momenta in order 
to respect the Lorentz invariance. In our case of nonrelativistic theory, on the other hand, 
corresponding symmetry is Galilean invariance and rotation symmetry. There is no obvious 
way of imposing a cutoff maintaining them. It is a rather annoying issue in applying the 
method to nonrelativistic cases. Fortunately, as we explain in the next section, the correct 
way of implementing a cutoff is clear in the two-nucleon system. In this section, however, 
we naively consider the flow equation defined in Minkowski spacetime and impose a cutoff 
only on three-momenta, though it breaks Galilean invariance. 

We begin by defining IR cutoff generating functional for nonrelativistic nucleons with 
source terms, i]^ and rj. The generating functional for this system is defined in terms of 
integration over Grassmann variables, N{x), by 

N^-R'^I^-N^I^^ {N\p))^ {R^a\p))^^ mp))^ , (3.2) 

where S'aq is an arbitrary bare action composed of local operators at physical UV cutoff 
scale Aq. The indices A and B denote spin and isospin collectively. Fourier transform is 
defined as follows. 

The function R^l\p) effectively cuts off the IR part of the fluctuations in the integrand so 



that M/a[?7''^, 77] contains the effects of only the fluctuations for momenta larger than A in 
the coupling constants. In principle, -Ra(p) can be an arbitrary function which satisfies the 
following properties, 

00 as (p/A)^ 0, 
as (p/A)2 00, 

but for definiteness we use the following IR cutoff function in our analysis^. 



(3.4) 



PL 

2M 









1 — exp 







(3.5) 



and take the sharp cutoff limit n ^ 00 in most cases. See the next section. 
For convenience we introduce the compact notation, 

Jn = {v,V^), Un = {N\N), (n=l,2), 
and write the generating functional 1)3.11) as 



(3.6) 



(3.7) 



where 



(-Ra^)ab 



(3.8) 



{R^^)ab 
{R^\)ab 

The Legendre transform rA[$] of 1Va[^] may be defined in the standard way by intro- 
ducing the expectation value of Un in the presence of the source J„, 

5Wk 



v^[^ = w^[j]-JA 



6J„. 



(3.9) 



where = has been introduced. 

It is more useful to define Fa, an averaged action, as follows, 

Ta = f A — -^n{RA)nm^Tn, 

which satisfies the following Legendre flow equation, 

dR\ 



(3.10) 



d£A ^ 
dA 2 ^ 



dA 



(r(2)+i?A)"' 



(3.11) 



^ In order to cutoff the momentum in an A/-independent way, we include M in the cutoff function in this 
form. 
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where we have introduced the notation {^{2))^^^ = S'^^A/^^n^^m and Tr denotes the inte- 
gration over momentum and also the trace in the internal space. Although Legendre flow 
equation resembles a one-loop equation, it contains nonperturbative information through a 
full propagator with IR cutoff, (r(2) + Ra) ^ 

We introduce a dimensionless parameter t = In (Aq/A) and dt which denotes the derivative 
with respect to t that acts only on explicit t-dependence of the IR cutoff function R^- The 
flow equation may be written in an even simpler form, 

^ = lTra[ln(r(,) + i?A)]. (3.12) 

We may further simplify the expression by using the nonrelativistic feature of the theory. 
For this purpose, let us split (r(2) + -Ra) into fleld independent part V and fleld dependent 
part JF. V is the "full propagator", while JF is the sum of the multipoint "vertices" which 
have two internal lines and $„ attached to each of the external legs. The RHS of ()3.12|) can 
be written as 



-Tr dt 
2 



(3.13) 



In nonrelativistic theory, there are no anti-particles so that the loop diagrams are very 
limited compared to relativistic theory; the diagrams with anti-particle lines do not exist. 
The flrst term in ()3.13p vanishes. (It is a fleld independent constant anyway.) The second 
term also vanishes because there is no way to draw the loop without an anti-particle line. 
The third and higher order terms contain the non-vanishing diagrams. 

In the following sections, we concentrate on the two-nucleon sector, to which only the 
four-nucleon (4N) operators contribute. It is easy to see that only the third term contains 
such contributions. (It also contains the contributions to the other sectors.) The other 
sectors get contributions from several terms in ()3.13p . To the three- nucleon sector, for 
example, both the third and the fourth terms contain the contributions. 

In this way, we end up with the following reduced equation for the 4N operators. 



= - -Tr dt [{V~':Fr] 

4N ^ 



(3.14) 

4N 



dt 

Although the above equation is exact, one needs an approximation to solve it. We consider 



a simple truncation of the space of operators. Namely, we consider only the operators of 
leading orders in derivative expansion. The approximation is based on our hope that, even 
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though some operators get large anomalous dimensions, their "ordering" of importance 
would not change; the lower the canonical dimension is, the lower the scaling dimension 
would be. Of course some mixing of operators would occur and should be taken into account 
properly. There is no absolute justification for this hope. One should examine its validity by 
actually enlarging the space of operators and confirming the stability of the results against 
the enlargement. 

An actual calculation goes as follows. We substitute the "ansatz" for the averaged action 
consisting of the operators up to a certain order into the above flow equation. The "ansatz" 
determines the explicit forms of V and JF. We then compare the coefficients of the operators 
on both sides. Higher oder operators which emerge in the RHS of (j3.14j) are disregarded 
("projected out"). The LHS contains the derivatives of the couplings while the RHS does 
not, thus the comparison leads to a set of first-order differential equations for the couplings, 
the KG equations. 

It is important to note that a WRG transformation generates all kinds of operators 
allowed by the symmetry of the theory including the so-called "redundant" operators, the 
operators which may be eliminated by the use of equations of motion. The "use of equations 
of motion" act ually means a field redefinition, which eliminates the "redundant" term. In the 
previous paperj27|, we showed that the field redefinition gives rise to a nontrivial Jacobian, 
and that "redundant" operators play an important role in a WRG analysis. 



IV. RG ANALYSIS OF PIONLESS NEFT 



A. pionless NEFT 

In this section, we consider the RG fiow of the pionless NEFT, in which only nonrela- 
tivistic nucleons are treated as explicit degrees of freedom. Effects of anti-nucleons, pions, 
and other heavier mesons and baryons are integrated out and hidden in the values of the 
coupling constants. Thus the theory is expected to describe the interactions of nonrela- 
tivistic nucleons with external momenta sufficiently smaller than pion mass. Although our 
approach is general and not restricted to this theory, it serves as the simplest example which 
illustrates the essential points. 

In a WRG analysis, symmetry is extremely important. We assume Galilean invariance, 
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rotational symmetry, and invariance under charge conjugation, parity, and time-reversal. 
Furthermore, we assume exact isospin symmetry, and ignore electromagnetic and weak in- 
teractions. In the following, we will focus only on two channels, ^5*0 and ^Si-^Di, in the 
two-nucleon sector. The extension to higher partial waves is easy. 

We study the RG flow for this theory by using Legendre flow equation. We consider 
the following "ansatz" for the averaged action, keeping only the operators with (canonical) 
dimension up to eight. 



' 2M 



[^Sq channel) 
[^Si-^Di channel) 



(4.1) 



where the operators in the ^5*0 are given by 



o 



(4.2a) 
(4.2b) 

(4.2c) 



and in the ^Si-^Di channel. 



(T) 



Of 



(4.3a) 
(4.3b) 

(4.3c) 

(4.3d) 
(4.3e) 



where we have introduced the notation = V^ + V^ — 2V-V and the projection 
operators [2^ ■ 

pi'so) = l^a^r^r\ pf = iaVV^, (4.4) 



for the ^5*0 channel and the ^Si channel respectively. The nucleon field N{x) with mass M 
is an isospin doublet nonrelativistic two-component spinor. Pauli matrices a* and r" act on 
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d 

dK 





four-nucleon vertices ® = ^^S^ 



dK 

FIG. 1: The loop diagram which contributes to the Legendre flow equation for four-nucleon oper- 
ators. 

spin indices and isospin indices respectively. The two channels are completely decoupled, 
and thus we can consider each channel separately. 

Note that we have not included the wave function renormalization because there is no 
such contribution due to the nonrelativistic feature. We neither include six-nucleon operators 
and higher because of the nucleon number conservation. The more-than-two-nucleon sectors 
are completely decoupled from the two-nucleon sector which we are interested in. 

Note also that we have included redundant operators, O2 and O2 , because they are 
necessary in a consistent WRG analysis. 

Let us consider the simplest example of the loop calculation in order to explain how 
we maintain the Galilean invariance. The diagrams which contribute to the Legendre flow 
equation are of the type shown in FIG.^ with various factors at the vertices. The propagator 
is IR cutoff by the function so that we are interested in the following integral, 

d'k ^ ^ ^^^^ 



(27r)4 po/2 + p - 7^ ((P/2 + ky) + te ((P/2 - kf) + le 

where we have introduced 

n {k') = ^ - {k^) , (4.6) 

and (P", P) is the momentum flowing in (and out of) the diagram. Note that we assign 
the loop momentum (/c", k) so that / is invariant under P —P. Note also that there 
may be additional momentum factors from the vertices but they do not affect the following 
argument. After integrating kP, we have 

<fk i 



(27r)3 po-n ((P/2 - kf) - n ((P/2 + kf) 



(4.7) 
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It is important to note that, if the IR cutoff were absent, the propagator would be 



^ 4M M ^ 



(4.8) 



where the combination P"^ /AM describes the center-of-mass kinetic energy of the two nu- 
cleons while Vt^ jM represents the kinetic energy of the relative motion. It is now clear that 
the only way to maintain the Galilean invariance is to cutoff the relative momentum. It is 
therefore physically legitimate to replace ()4.7|) with 

[ J-L I (4 9) 

J (2.)3po_^_g + 24^)(fc^)+.e- ^ • ^ 

The integral may be evaluated with an arbitrary parameter n. A useful formula for a 

typical loop integral is given in Appendix lA II In the following, we consider only the sharp 

cutoff limit n ^ oo, for which the RG equation becomes the simplest. We also look at the 

dependence of the fixed points and the anomalous dimensions against the variation of n. 

They are shown in Appendix lA 21 



B. RG flows and flxed points 



It is useful to define the dimensionless variables as follows and write the RG equation in 
terms of these variables. In the ^5*0 channel, we introduce^, 

27r2 



^^^cr, y 



27r2 



2 5 



27r2 



in terms of which the RG equations are given by 



dx 
'dt 



—X — 



dy o 

-dt = -'y- 

dz 



x^ + 2xy + y"^ + 2xz + 2yz + z^ 



^x^ + 2xy + ^y^ + yz - ^z^ 



^x^ + xy + ^y^ -xz-yz- ^z^ 



(4.10) 



(4.11a) 
(4.11b) 
(4.11c) 



^ The factor M comes from the scale transformation property in nonrelativistic theory, t' = X^t, x' = 
Xx, where A is the scale factor. Under this transformation, the nucleon field transforms as N{t,x) = 
X^ N' [t' ,x'), and thus we have 00^"* = X^O'\^^ for example. The d'^x gives the factor A~^, and the 
additional factor is supplied by the A-dependence of the coupling Cg'^'' oc to cancel the A^. To 
have the right mass dimension, Cq has the 1/M dependence. Similar consideration leads to the correct 
A and A/ dependence of other couplings. 
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Note that the RG equations are quadratic in four-nucleon couphngs due to the nonexistence 
of anti-nucleons in our nonrelativistic formulation, up to the first terms, which come from 
the canonical scaling. 

Similarly in the ^Si-'^Di channel, we introduce 



27r2 ° ' ^ 27r2 
and we have the following RG equations 



27r2" 



27r2 3 



"2 ! 



(4.12) 



dx' 
'dt 



-X 



dy' o / 

-di = -'y 



dz' 
dt 



dw' 
~dt 



-?>w' 



x'"^ + 2xy' + y'"^ + 2x z + 2yz + z'^ + 2w'^ 



1 /2 , ^ / / 3 ,2 , / / 1 /2 , /2 

-X + 2x y + -y +y z - -z +w 



1 /2 , / / 1 /2 / / / / 3 ,2 , ,2 

-X + X y + -y — xz— yz— -z +w 



II, II, II 

X w + y w + z w 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 



Note that the RG equations ()4.13|) for the ^Si-^Di channel are identical to the RG equations 
(14. lip for the ^5*0 channel if w' is set equal to zero. (This is actually a solution.) 

We can now draw the RG flow diagram. FIG. El shows the RG flow for the ^5*0 channel 
and FIG. 01 for the ^Si-^Di channel up to operators with dimension eight. Both flows are 
projected on to the two-dimensional plane spanned by the lowest order coupling constants. 
As expected from the RG equations ()4.11|) and ()4.1 



the RG flows for the ^Sn channel and 



for the ^Si-^Di channel have similar structure. First of all, there are two fixed points in 
both RG flows; the trivial IR fixed point and the non-trivial UV fixed point. Secondly, both 
RG flows have two phases; one in which all points approaches the trivial fixed point in the 
limit of vanishing cutoff, and the other one in which all points approaches infinity at certain 
finite cutoff scale. The nontrivial fixed points are on the phase boundaries. 

Location of the fixed points are given by solutions of the equations, dx/dt = dy/dt = 
dz/dt = 0, for the ^5*0 channel, and dx' /dt = dy' /dt = dz' /dt = dw' /dt = for the ^Si-^Di 
channel. As explained in Appendix IA21 we find four fixed points in each channel for an 
arbitrary value of the parameter n. One of them is extremely unstable against the variation 
of n, and may be considered as "spurious", i.e., it is an artifact of our truncation. The rest 
are stable and take the following values in the n —>■ oo limit. 
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FIG. 2: The RG flow for the ^Sq channel 

(S) (S) 

projected onto the Cq -C2 plane. Flow- 
lines start with the arbitrarily chosen value 
2; = at the edges of the graph. 

^5*0 channel: 



FIG. 3: The RG flow for the ^Si-^D^ 

(T) (T) 

channel projected onto the Cq -C2 plane. 
Flow lines start with the arbitrarily chosen 
values z' = and w' = 0.1 at the edges of 
the graph. 



{x\f,zn = (0,0,0), (-1,-^,0, (-^'f'-^) ^^-^^^ 

^Si-^Di channel: 

y'\ z'\ w'*) = (0, 0, 0, 0) , (^-1, ~, ^, 0^ , (^-9, y , 0^ (4.15) 

We also disregard the third fixed points in both channels, because they are irrelevant to 
the discussion of the power counting. The eigenoperators at these fixed points are complex 
and have complex eigenvalues (scaling dimensions). It means that the operators do not have 
the definite scaling property there. 

We obtain the scattering amplitudes explicitly in both channels in Appendix |Bl as we 
did in the previous paper j22|, and find that the RG equations, which look much more 
complicated than 1)4.111) and ()4.13|) . have only two fixed points in each channel and they 
exactly agree with the first two fixed points obtained here'*. 

^ Actually, the RG equation obtained in the previous paper^ is those for the ^5*0 channel. 
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C. scaling dimensions 



We now study the behavior of the couphng constants around each fixed point by consid- 
ering the hnearized RG equations. 

Let us first consider the trivial fixed points. In the ^Sq channel, the trivial fixed point is 
at {x*,y*,z*) = (0,0,0). We consider a small deviation {6x,6y,6z) around it and linearize 
the RG equations ()4.1H) . 



d 
dt 



5y 



( 



\ 










-3 
-3 



by 
5z 



(4.16) 



The eigenvalues and corresponding eigenvectors are trivially obtained, 

1^2 = -3 : U2 = 



-1 : Ui 





VV 



-3 : 





VI 



(4.17) 



Namely, the couplings C'^\ Cf\ and B^^^ have the scaling dimension —1, —3, and —3 
respectively. This is what one naively expects from their canonical scaling properties. 

Similarly in the ^Si-^Di channel, we consider a deviation (5x', 5?/', 5z' ^ 6w') around the 



trivial fixed point (x' 
tions from Eq. (j4.13|) 



y ,z ,w ) 



(0, 0, 0, 0), and find the following linearized RG equa- 





/ 6x' \ 


1 


-1 








\ 


/ 5x' \ 


d 


5y' 







-3 








by' 


dt 


5z' 










-3 





6z' 




\6w' J 













-3/ 


\6w' J 



(4.18) 
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with the eigenvalues and corresponding eigenvectors, 



-1 : Ml 



Z/4 



-3 : Ui 






/o\ 




VI/ 



z/2 = -3 : U2 



1 




= -3 : Ms 





1 

vo/ 



(4.19) 



The couphng constants C'^\ B^'^\ and €2^^^ have the scahng dimension —1, —3, —3, 
and —3 respectively. 

Note that all of the scaling dimensions are negative at the trivial fixed point and thus 
the corresponding operators are irrelevant, i.e., they should be treated perturbatively. 

At the nontrivial fixed point, on the other hand, something much more interesting hap- 
pens. In the ^5*0 channel, the deviations from the nontrivial fixed point {x*,y*,z*) = 
(— 1, —1/2, 1/2) satisfy the following linearized RG equations. 



d 
dt 



y6zj 



/ 1 2 2 
2 1 
y-2 -2 -3 



5y 



(4.20) 



with the eigenvalues and corresponding eigenvectors. 



ui = +1 : ui 



( ' ] 

1 



, U2 = -1 : U2 



I ^ 

-1 

V 1 ) 



z/3 = -2 : U3 



/ 2 N 



(4.21) 



Note that the scaling dimensions change drastically. Note also that the eigenoperators are 
not monomials but linear combinations of the operators. This is a typical feature of the 
momentum cutoff regularization. 

Similar results are obtained in the ^Si-^Di channel, where the linearized equations around 
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the nontrivial fixed point {x'* , y'* , z'* , w'*) 

/ 6x'\ / 1 



d 
It 



5y' 
5z' 
\6w' J 



(-1,-1/2,1/2,0) become 
2 2 \ / 5x'\ 



2 10 
-2 -2 -3 
\ -2 / 



6y' 
6z' 
\6w' J 



(4.22) 



giving the following eigenvalues and corresponding eigenvectors. 



1^1 



1 : Ui 



1^3 



-2: U3 



( 2 \ 

-1 
-2 
V / 



Z/4 



M4 



(4.23) 



/ 1 \ / \ 

1 -1 

, Z/2 = -1 : U2 = 

-1 1 

V / V / 

/o\ 





VI/ 

Note that the structure is essentially the same as that in the ^6*0 channel, with bw' decoupled. 

The most important result at the nontrivial fixed point is that there is a relevant coupling 
in each channel. The coupling ui has dimension +1, and according to the argument given 
in Sec. ini the corresponding operator should be treated nonperturbatively, while the others 
perturbatively. The operators should be ordered according to the scaling dimensions. The 
less the dimension of a coupling constant is, the less the importance of the corresponding 
operator is. 

As we have done in the previous paper[27|, RG equations may be obtained by directly 
calculating the amplitudes and requiring them to be independent of the "floating" scale A. 
See AppendixElfor the results. The resulting RG equations are much more complicated. It is 
interesting, however, the location of the fixed points and the scaling dimensions are precisely 
the same as those obtained here with the Legendre flow equation. The eigenvectors are a 
bit different, though. 



D. physical meaning of the RG phases 

Let us consider the real world and examine the physical meaning of what we have found 
in the two-nucleon system. 
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As mentioned in Introduction, the two-nucleon system is considered as "fine-tuned," in 
the sense that the scattering length is unnaturally large. It suggests that the actual system 
is near the nontrivial fixed point in each channel. The closer to the fixed point it is, the 
longer the scattering length is. 

What then makes the physical difference between the ^5*0 and the ^Si-^Di channels, both 
of which are close to the nontrivial fixed points? The answer to this question is provided 
by the fiow diagram; the nontrivial fixed point is on the boundary between two phases. In 
one of them, the fiow goes toward the trivial fixed point in the IR. We therefore call this the 
weak coupling phase. In the other, the fiow goes toward stronger couplings, thus we call it 
the strong coupling phase. It is natural to think that the real two-nucleon system in the ^Sq 
channel, in which there is no bound state, is in the weak coupling phase, while that in the 
^Si-^Di channel, in which there is a shallow bound state, deuteron, is in the strong coupling 
phase. 

In order to make the physical picture more transparent, it is useful to examine the four- 
nucleon (two-nucleon scattering) amplitudes. As shown in Ref. and in Appendix El 
the four-nucleon amplitudes may be obtained explicitly in the derivative expansion. By 
substituting the solution of the linearized RG equation in the ^5*0 channel obtained from the 
direct method, 



y6zj 



( 2 \ 



v-v 



Ao 



/ ^ 



V 1 1 



A 



( 1 \ 

1 

v-v 



Ao 
A 



(4.24) 



where a, 6, and c are small dimensionless constants, into the amplitude in the center-of-mass 
frame, we obtained the renormalized (off-shell) amplitude near the nontrivial fixed point [2^, 



M 
47r 



A^\p\Ti\p' 

where ellipsis denotes higher orders in I/Aq. 

By comparing it with the effective range expansion 

1 1 



TT vr V Ao ' 



+ 



on— shell 



M 



_l — fp — ip 
a 2 



(4.25) 



(4.26) 



on the mass shell p = Mp^ = \p\ = one sees the scattering length a and the effective 



21 



range r are given by 

a = -^, r = -^. (4.27) 
2cAo ttAo ^ ' 

That is, the inverse of the scattering length is identified as the relevant coupling. The fine- 
tuning |a| ^ Aq^ indeed corresponds to |c| ^ 1. It is also important that the sign of the 
couphng c distinguishes the phases, namely, the inverse of the scattering length is the order 
parameter. It vanishes on the boundary, the critical surface. Because the effective range is 
of natural size, r ~ Aq ^, it is more accurate to say that the real two-nucleon system is close 
to the critical surface. 

The location of the pole which is found within the range of EFT is given by 



A2 

^^0 



(?)■ 


166c" 







(4.28) 



A similar expression is obtained for the pole in the ^Si-^Di channel. The pole on the 
physical Riemann sheet that corresponds to a bound state occurs only when c < 0, i.e., in 
the strong coupling phase, while there is no restriction on b. 



V. COMMENTS ON OTHER POWER COUNTING SCHEMES 

As we emphasized in Sec. m our way of determining the power counting on the basis of 
the scaling dimensions is general and nonperturbative. It does not depend on any specific 
regularization/renormalization scheme, nor on which fixed point we are looking at. 

In this section, we would like to comment on other works from our point of view and make 
the connection between their power counting and ours. It leads us to a deeper understanding 
of the power counting issues. 

First of all, we would like to mention the work done by Birse, McGovern, and 
Richardson [2^, which is very close to ours in philosophy. They consider the most gen- 
eral potential and the Wilsonian RG equation satisfied by it. The "fixed point" potential 
depends only on the energy, 

(5.1) 

where p = V ME, with the E being the center-of-mass energy. They then consider the 
perturbation around the "fixed point" potential, and find the eigenvalues of the linearized 
RG equations. The first few eigenvalues coincide with ours. They use the information to 
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Vo{p,A) 



A-^ln^ 
2 A-p 



determine the power counting. It is important to note that the "energy-dependent perturba- 
tions" have the spectrum u = —1,1,3, ■■■ , while the "momentum-dependent perturbations" 
have u = 2,4,5, ■■■ . (Note that their u^s have opposite signs to the scahng dimensions used 
in this paper.) 

To compare their resuhs with ours, it is useful to translate our operator formulation to 
the potential one. Our operators may be equivalent to the potential, 

.2 I _/2- 



Vip',p,p') =Co + AC2ip'+p") - 25 - + 



(5.2) 



where we suppress the superscripts which denote the channel. If we substitute the fixed 
point values of the couplings, we have 



V*ip',p,p') 



271 



2 



1 + 



(5.3) 



MA 

which is exactly the expansion of the potential ()5.1|) . with p^ replaced by E. Note that the 
momentum dependence cancels at the fixed point. This may be an evidence that both ap- 
proaches are essentially the same. Actually, the "perturbations" with the scaling dimensions 
z/i = 1 and z/2 = — 1 corresponds to the potentials that depends only on while that with 
1/3 = —2 to a momentum-dependent one, in agreement with their findings^. 

For this simplest pionless NEFT in the two-nucleon sector, the method based on the 
potential has an advantage. On the other hand, because our approach is based on the general 
principle of fully-fledged quantum field theory, we could treat more complicated systems in 
a unified and systematic way. For example, we could treat three-nucleon problems without 
worrying about the off-shell ambiguities, and fully relativistic systems for which the potential 
picture is not valid. The important point is that our work gives a field-theoretical foundation 
for their approach and has opened a window to wider applications. 

In the following, we will make comments on other power counting schemes on the basis 
of what we have found. 

In Weinberg's schemey,^!, one assumes the NDA power counting in making the "effective 
potential," then treats it nonperturbatively. As we have explained, however, the NDA is the 
power counting near the trivial fixed point, where all the operators are irrelevant and should 



^ A careful study reveals that the eigenvector for 1/3 = —2 does not agree with theirs, though those for 
1/1 = 1 and 1^2 = —1 do. The reason seems that the reaction matrix they consider does not have a direct 
connection to the scattering amplitude which we work with off the mass shell. 
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be treated perturbatively. The so-called infrared enhancement due to two-nucleon reducible 
diagrams somehow remedies the above-mentioned mismatch in the treatment; it makes the 
"effective potential" more relevant. The way is not systematic however, and in some cases 
inconsistent . 

The scheme due to Kaplan, Savage, and WiseQ, Q is important because they em- 
phasized one of the most important aspects of the two-nucleon system; it is close to the 
nontrivial fixed point! (The existence of the nontrivial fixed point itself had been noticed 

ri 

by Weinbergp].) It is however very subtle why the KSW scheme succeeds in some channels, 
and fails in others. To understand the points, we have to examine it very carefully. 

It is difficult to define the KSW power counting without mentioning the Power Divergence 
Subtraction (PDS) scheme®. In PDS, the poles at Z) = 3 as well as -D = 4 are subtracted. 
In the pionless NEFT, there is no pole at D = 4, so that one subtracts as if one were 

in D = 3 dimensions. This is the heart of the PDS scheme; it changes the (canonical) 

( s) 

dimensions of the operators! For example, the operator Oq has the canonical dimension 
6 in D = 4 dimensions, but would have 4 in = 3 dimensions. In PDS, it is treated as 
an operator of dimension 4. This shift of the dimension by two coincides with the correct 
anomalous dimension of the operator, as we found in the previous section. (Remember that 
the leading contact operators have anomalous dimension 2, and according to Birse et al, all 
the "energy-dependent perturbations" seem to have the same anomalous dimensions.) We 
believe that it is the true reason why the KSW/PDS scheme works near the nontrivial fixed 
point^. 

Actually, this "shift by two" applies to all the four-nucleon operators in the KSW/PDS 
scheme^. To see this, consider the KSW/PDS RG equation for the dimensionless couplings. 



where 72^ = (M/i2"+747r)C2„, with C2„(/i) being given in Eq. (2.10) in Ref. (l9|. Near the 



^ We do not claim that PDS is indispensable for the KSW power counting, but that the PDS scheme 

provides the most unambiguous definition of the KSW power counting. 
^ From this point of view, it is legitimate to think that the KSW/PDS scheme should be applied only to 

the case near the nontrivial fixed point, though the RG equations (|5.4|l know about the trivial fixed point 

too. 

® The redundant operators scales differently than the others, but they are not necessary with dimensional 
regularization. 




(5.4) 
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nontrivial fixed point (7q, 72, 74, ■ ■ ■ ) = (—1, 0, 0, ■ ■ ■ ), the RHS of ()5.4p becomes 

(2n + 1)72„ + 27o72n + ■ ■ ■ ~ (2n - 1)72^, (5.5) 

wliere tlie ellipsis stands for terms quadratic in 7's and independent of 70, which do not 
contribute to the anomalous dimensions. This explicitly shows the shift from {2n + 1) to 

As far as all the anomalous dimensions are two, the KSW/PDS scheme would work 
perfectly. The problem in the pionful NEFT seems that the singular interactions due to pion 
exchange in the ^Si-^Di channel would modify the anomalous dimensions of the contact 



operators strongly, making the "shift by two" rule invalid. Actually, Birse[2j] found that the 
anomalous dimensions shift additionally by minus one half at the nontrivial fixed point in 
the presence of the tensor force due to one-pion exchange. The investigation in this direction 
is currently in progress [3l|. 

Several authors considered power countings which are compatible with the KSW power 
counting without using PDS so that one might think that the KSW scheme may be defined 
independent of PDS. In fact, it appears possible to do so if one assumes the scaling of the 
coup ling constants which is a consequence of the PDS regularization/renormalization. In 
Ref. j32|, van Kolck introduced a scale K, which accounts for the fine-tuning, and assumed 
the scaling of the coupling constants without specifying the regularization/renormalization. 
The scale K seems to correspond to our cAq with c -C 1 introduced in Eq. ()4.24|) . It is 
important to note that K is a physical scale so that there is no notion of RG with respect to 
K in his approach. If one nevertheless identifies N with the scale fi in the KSW scheme, the 
assumed scaling of the couplings is just the one in the KSW/PDS scheme. On the other hand, 



using momentum space subtraction renormalization, Mehen and Stewart 
renormalization condition (Eq. (6) in Ref. 
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imposed the 



^^^'""'^ l™,..^o = -^C,M i^l^nr , (5.6) 



^ There is a subtle point about the scaUng of the couphng constants. If one solves the RG equations near the 
fixed point, the couphng constant scales C2n ~ ^tt/MAq"'^^"^, but it is different from the KSW scahng 
C2n ~ 4:TT / M Aq fi^~^^ . It is due to the assumption that the effective range is of natural size, r ^ I/Aq. 
Namely, the KSW scaling is not that for a point close to the fixed point, but for a point close to the 
critical surface. 
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on the (part of the) amphtude that scales as Q"^~^, where Q is a typical momentum 
scale. (This power counting for the amplitudes is also assumed in accordance with the 
KSW/PDS scheme.) Putting /i/j ~ Q, it requires the coupling constant C2m{l^R) to scale as 
~ Q^~^ (/^j?)"^"* ~ Q~^'"^~^^\ which is essentially equivalent to the KSW/PDS power count- 
ing. It is still unclear to us why the renormalization condition (|5.6p reflects the situation 
that we are close to the nontrivial fixed point. 

Let us finally comment on the confusion in the literature over the scale introduced in 
dimensional regularization. It is important to note that the scale does not have intrinsic 
meaning of "typical scale." It is completely arbitrary and physics is independent of it. In 
the usual applications of dimensional regularization to relativistic field theory, /i appears in 
logarithms as ln((5^//i^), where Q is the typical scale of the process. It is this dependence 
that makes it useful to think /i ~ Q in order to suppress the "large logarithms." In our 
application to NEFT, however, there is no logarithm, thus no reason to think /i ~ Q. 

VI. CONCLUSION 

In this paper we studied the Wilsonian RG flow of pionless NEFT using the Legendre 
flow equation in order to determine the correct power counting. We emphasized that the 
determination of power counting should be based on the scaling dimensions, because power 
counting is in essence the order of magnitude estimate based on dimensional analysis, and 
the scaling dimensions are quantum dimensions of the operators. 

The RG flow has a nontrivial fixed point as well as a trivial one in each of the ^Sq and 
^Si-^Di channels. Near the nontrivial fixed point, we found a relevant operator, the coupling 
constant of which is proportional to the inverse of the scattering length. The nontrivial fixed 
point is on the boundary of two phases, the weak coupling and the strong coupling phases. 
The inverse of the scattering length is regarded as the order parameter, which takes zero 
at the boundary. The real two-nucleon systems are considered to be close to the nontrivial 
fixed point. 

The difference between the ^5*0 and ^Si-^Di channels comes from the difference of the 
phases they are in. In the strong coupling phase, the scattering length is positive and the 
four-nucleon amplitude has a physical pole which represent a bound state, the deuteron in 
the ^Si-^Di channel, while in the weak coupling phase, the scattering length is negative and 
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there is no bound state; this phase describes the ^5*0 channeL We beheve that this way of 
viewing the difference between the two channels is new, and we expect that it would give 
some further insight into the character of nuclear force. 

The relation to the work by Birse et al. was clarified. The use of fully developed EFT 
framework was emphasized. We also made comments on other power counting schemes and 
tried to understand why they succeed in some cases and fail in others from our point of 
view. In particular, we claimed that the true reason why the KSW/PDS scheme works is 
identified as the "shift by two," namely, PDS happens to capture the essential feature that 
the anomalous dimensions of the contact two-nucleon interactions are two. 

It was essential to include the so-called redundant operators to consistently perform the 



Wilsonian RG analysis 



27| in a fully field theoretical manner. They are also necessary to 



renormalize the off-shell amplitudes. The inclusion of them seems essential to consider the 
tree-nucleon systems. 

The RG equations derived from the Legendre flow equation are much simpler than those 
from the amplitude, while the important physical information (the existence and the location 
of fixed points, the global structure of the flows, the scaling dimensions, etc.) is the same. 
The use of the Legendre flow equation would help us to tackle more difficult problems, where 
the direct method seems infeasible. 

It is very important to apply the present method based on the Wilsonian RG to determine 
the power counting of the pionful NEFT, where the power counting issues are most contro- 
versial. The singular tensor force would change the anomalous dimension of the lowest-order 
contact operator drastically in the ^Si-^Di channel, invalidating the "shift by two" rule. One 
might imagine that, at a high energy scale, there is actually no difference of the values of 
couplings between the ^5*0 and ^Si-^Di channels, but the pion interactions make them so 
different. This scenario appears attractive from the approximate spin- flavor symmetry point 
of view Note that the nontrivial fixed point is spin-flavor symmetric. 

One of the interesting results in this paper is the explicit relation between the Wilsonian 
RG analysis and the existence of the bound state in the strong coupling phase. It would 
open a new possibility of studying bound states using the Wilsonian RG. On the other hand, 
it would become important when we study the tree-nucleon systems without introducing a 



"dibaryon" field jS^ 
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40|. Note that most of the existing literature on the NEFT study 



of three-nucleon systems is based on the "dibaryon" formalism. Because our formulation is 
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fully field theoretical, it is free from the off-shell ambiguities which the approaches based on 
potentials suffer from. 



Acknowledgments 

We would like to thank Y. Yamamoto for several useful discussions. We are also grateful 
to H. Gies for pointing out the wrong placing of the derivative dt in Sec. 11111 in the earlier 
manuscript. One of the authors (K.H.) is partially supported by Grant-in-Aid for Scientific 
Research on Priority Area, Number of Area 763, "Dynamics of Strings and Fields," from 
the Ministry of Education, Culture, Sports, Science and Technology, Japan. 



APPENDIX A: IR CUTOFF FUNCTION WITH FINITE n 

In this Appendix, we give a useful formula involving the IR cutoff function, and the 
dependence of the fixed points and the scaling dimensions of the RG equations on n. 



1. Evaluation of a typical loop integral 

A typical loop integral with the IR cutoff function R^^\k'^) which appears repeatedly in 
the derivation of our RG equations is 

(2vr)3^*^(p)_g + 2i?«(fc2)' ^ ^ 

where A[P) is a Galilean invariant quantity, 

and = (P°, P) is the momentum flowing in (and out of) the diagram, P = pi + P2 = 
P3+P4, where pi and p2 are initial momenta, and ps and p^ are final momenta of the nucleons. 
The IR cutoff function is chosen as (j3.5p for definiteness. 

Switching to the polar coordinate, making the change of variable to a; = k/A, and ex- 
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panding the integrand with respect to A{P) = (M/A^) A{P), we have 



^ Jo 

2m + 1 \ ~, , I 
A{P)T 1 



R{x) 



MA 



2m+l 



27r2 



r 1 



2n 



A{P) -x^(l- R{x)^ 

2m — l\ / l-2m 

2 _ 



2/2 



(A3) 



where R{x) = (1 — exp (x^")) ^ has been introduced. In the present paper, we disregard the 
terms of O^A"^) as higher orders. 

We may decompose A{P) into the following Galilean invariant quantities, 



A{P) 



o_ Pi 
2M 



P2 

P2 



2M AM 



^ (Pi - P2f 



(A4) 



Combining possible multiplicative momentum factors outside the integral, we can identify 
the contributions to the various operators in our theory space. For example, contributions 
from the first two terms renormalize a redundant operator and the last term does the oper- 
ators which contain spatial derivatives. 



2. Dependence on the IR cutoff function 

The RG equations are derived and solved with a finite value of n. In this section, we 
show the n dependence of the fixed points and the scaling dimensions. 

The dependence on n might give some useful information on the behavior of the present 
approximation. Note that if no approximation were made, the universal quantities such as 
the scaling dimension would not depend on how the RG transformation is defined. Weak 
dependence would suggest that the approximation is good. We will see shortly that the 
most of the fixed points as well as the scaling dimensions have very weak dependence. The 
fixed point with strong dependence on n may be regarded as "spurious" and, for other fixed 
points, the approximation is very good. 

In general, a sharp momentum cutoff gives a poor convergence and a smooth cutoff 
function is preferred j4l|. In the present case, however, the dependence on n is very weak 
and the use of a sharp cutoff is expected to be reasonable. In addition, it gives the simplest 
expressions. 
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TABLE I: The dependence on n of the fixed point A. 



n location of the fixed point(rE*, y*, z*) scaling dimensions 

2 (-1.10326,-0.60467,0.60467) (-2,-1,1) 

10 (-1.02722,-0.524991,0.524991) (-2,-1,1) 

10^ (-1.00287,-0.502587,0.502587) (-2,-1,1) 

10^ (-1.00029,-0.50026,0.50026) (-2,-1,1) 

10^ (-1.00003,-0.599926,0.500026) (-2,-1,1) 

oo (-2,-1,1) 



For a finite value of n, we find four fixed points including the trivial one. The dependence 
of the three nontrivial fixed points in the ^5*0 channel is given in Tables HI ITTl and IIIII (The 
dependence is similar for the '^Si-^Di channel.) The last column shows the three eigenvalues 
of the linearized RG equations at the fixed point. Let us call these fixed points A, B, and 
C. 

Note that, although we have an analytic expression for the RG equations for any n, 
the actual evaluation is done with the aid of a computer. We use Mathematica, but some 
numerical errors are not avoidable. 

One can easily see that the fixed points A and B have little n dependence (for n > 10), 
while the fixed point C depends on n strongly. We disregard the fixed point C. Note that 
even for the smallest integer value of n, n = 2 (n = 1 does not satisfy the condition ()3.4|) .). 
the scaling dimensions for the fixed point A are the same while the fixed points B and C 
disappear there. 

The fixed point B has complex eigenvalues, though they have weak dependence on n. As 
discussed in the main text, the eigen-operators do not have the definite scaling properties 
there and we disregard it too in our discussion of the power counting. At the moment we 
know neither if it is just an artifact, nor, even if it is not, its physical implications. 
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TABLE II: The dependence on n of the fixed point B. 


ft 


loration of the fixed DoiTiti t* ?/* 


scalinp" diiTiensioris 


z 


Does not exist 


Does not exist 


10 


(-9.85064, 8.33328, -1.27074) 


(-0.79483 ±2.18601i, 3) 


102 


(-9.03165, 7.52438, -1.46746) 


(-0.502903 ±2.396H, 3) 


10^ 


(-9.00266, 7.50188, -1.49663) 


(-0.50003 ±2.3979i, 3) 


10^ 


(-9.00026, 7.50018, -1.49966) 


(-0.494178 ± 2.4005i, 3.00465) 


oo 


f_9 15 _3^ 


(-i±z^,3) 


TABLE III: The dependence on n of the fixed point C . 


n 


location of the fixed point(x*, y*, z*) 


scahng dimensions 


2 


Does not exist 


Does not exist 


10 


(-71.9608, 77.804, -5.93769) 


(-62.929,0.874861,3) 


10^ 


(-6226.88, 6961.68, -702.804) 


(-6202.06,0.99871,3) 


10^ 


(-609432, 681363, -71583.4) 


(-609248, 0.999984, 3.81652) 


10^ 


(-6.08078 X 10^6.79826 x 10^-7.17129 x 10^) 


(-6.0806 X 10^1,3) 


oo 


Does not exist 


Does not exist 



APPENDIX B: EXPLICIT DERIVATION OF RG EQUATIONS FROM THE AM- 
PLITUDES 

In this Appendix, we consider the (off-shell) two-nucleon amplitudes in the '^Si-^Di 
channel in order to derive the RG equations. (The amplitude in the ^5*0 channel is given in 
Ref. j^lO The amplitudes take a matrix form due to the presence of the ^Si-^Di mixing 
and satisfy the following Lippmann-Schwinger equation in the center-of-mass frame, 

-^A^j{po,pl,pfj = -iVij {po,pI,pI) 

/d^k i 
jTT-r^ {-iVik {po,pI fc')) p {-tAkj {po, k^,pl)) , 

(Bl) 

where the subscript labels the partial waves, i,j,k = {S,D}, and Vij stands for the vertex 
in momentum space. The momentum integral is taken over < A; < A. 
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We take the averaged action ()4.ip as our action, and write the vertex Vij to order 0{p'^ 



Vss = cP + AC^ (pl + pl) - 2B(-) - 



SD 



(SD) 2 _ 4V^^(5D)_2 

"<-^2 Pn ^^ilS — — 



'^DD = 0, 



(B2) 

(B3) 



3 " ^ 3 
where is the center-of-mass energy, pi and p2 are the relative momenta in the initial 
and final states respectively. We here suppress the spin and isospin dependence. See the 
appendix in Ref. j20| for a useful technique. 
We consider the following ansatz. 



^55(/,P2,P?) = n{po)+T2{po){pl+pl)+Ts{po)plpl 

Asd{p\pIpI) = tMpI + tmpIpI 

ADsiPo,pl,pl) = Tg{po)pI + T7{po)pIpI, 

Add{po,pI,pI) = T8{po)pIpI, 

and obtain the solution, 



T3 



1 

D 
1 

D 



-(4a 



(T) 



M 



4C 



(T) 



M / 



M 

1 4v^^(5z^) 



D 3 



^7 = 

32 1 



1 4v^^(5D) 



D 3 



1 - 

AC. 



(T) 



M J 



Vo, 



a 



iSD)\' 



(B4a) 
(B4b) 
(B4c) 
(B4d) 

(B5a) 
(B5b) 
(B5c) 
(B5d) 

(B5e) 
(B5f) 



with 



D = l- (Cf) - 2Bi-V) /o - 2 (acP + ^) /i + (4Cf + ^) ' 



- I ( 4Cf ) + 



where are integrals defined by 

M 
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iSD)\ 



hi 



0-'2, 



27r2 



A /^2n+2 

dk 
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/i 



^/-Mp^ 



le. 



(B6) 



(B7) 



The RG equations for the couphngs are obtained by requiring that the inverse of the off- 
shell amplitudes {Ass,^sd in this case) are independent of the cutoff A, in the expansion 
in powers of pf, and /i^. We introduce dimensionless coupling constants as (|4.12p . 
The resulting RG equations are written as 



dY 



Y 



A — = — (QX^Z -3X^Z + 2XYZ + Y'^) 



A 



dA 
dw' 

llA 



X2 
1 

X 



{3X^w' 



Yw') 



where we have introduced X = 1 + {y' + z') /3, Y = x' — ((y' + z' 
2y' + ((y' + 2')^ — 2w'^) /3. Compare them with those obtained in Ref. 
channel. They are identical if one sets w' = 0. 



(B8a) 
(B8b) 
(B8c) 
(BSd) 

^ + 2w'£L/5, and Z 



for the ^5*0 



Although they look complicated and very different from Eqs. ()4.13|) . when expanded in 
couplings, they agree with Eqs. ()4. 13|1 to the second order. 

These RG equations have only two fixed points (as far as Mathematica can find); a trivial 
fixed point {X*,Y*, Z*,w*) = (1,0,0,0) and a nontrivial fixed point {X*,Y*, Z*,w*) = 
(1,-1,-1,0) which corresponds to {x'* , y'* , z'* , w'*) = (-1,-1/2,1/2,0). The locations of 
these fixed points are the same as those obtained from Eqs. ()4.13|) . 

Substituting 



x' = x'* + Sx', y' = y'* + Sy', z' = z'* + 6z', w' = w'* + 5w', 
into ()B8|1 . we obtain the linearized RG equations at the nontrivial fixed point. 



d_ 

dt 



(B9) 



/ 6x' \ 


/ 


1 


2 


2 


\ 


/ 6x' \ 


6y' 




2 


2 
3 


5 
3 





6y' 


Sz' 




-2 


8 
3 


11 
3 





5z' 


\6w' J 


\ 











-2/ 


\6w' / 



(BIO) 



where we have introduced t = \n (Aq/A). 
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These equations can be solved easily. 

f 5x'\ / 2 \ 

6y' ^ . 1 
5z' 
\6w' J 



+ b' 



4 

V / 

/ \ 

-1 
1 

V / 



A 



A 



/o\ 




VI/ 
/ 1 \ 

1 

-1 
V / 



A 

a; 



Ao 
A 



(Bll) 



M 
An 



+ ■ 



(B12) 



where a[, a'21 and d are dimensionless constants. This corresponds to Eq. (j4.24p . 

We obtain the renormalized off-shell amplitudes near the nontrivial fixed point to the 
second order in powers of and by substituting the solutions of linearized RG equations, 
for example, 

"2 4 /A/po~ 

-cAo - -6 — 

TT TT V Ac 

where the ellipsis denotes terms of higher orders of P^/Aq and p^/Ap. We notice that a'^ and 
a'2 do not contribute to the inverse of the amplitude to this order. The existence of such 
eigenvectors whose coefficients do not show up in the on-shell amplitudes is related to the 
equivalence theorem considered in the previous paper j23]. Actually, the equivalence relation 
(3.21) of Ref. [s^] gives "trajectories" in the space of couplings, and the tangential direction 
at the nontrivial fixed point is the direction of the eigenvector whose coefficient does not 
appear in the physical amplitude. This would suggest that all the "momentum-dependent 
perturbations" discussed in Ref. j^l form "surfaces" in the space of couplings near the fixed 
point on which any two points are related by some equivalence relation. 
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